#include<iostream>
#include<vector>
#include<cstring>
using namespace std;
const int N=5e2+10;
int n,m,s,d;
int dis[N],a[N],g[N][N],vis[N];
vector<int>res,per;
int sum,cnt;
void bfs()
{
	int res=0;
    for(int i=0;i<n;i++)
        dis[i]=g[s][i];
	dis[s]=0;
	for(int i=0;i<n;i++)
	{
		int cur=-1;
		for(int j=0;j<n;j++)
		{
			if(!vis[j]&&(cur==-1||dis[cur]>dis[j]))
			cur=j;
		}
		vis[cur]=1;
		for(int j=0;j<n;j++)
		{
			dis[j]=min(dis[j],dis[cur]+g[cur][j]);
		}
	}
}
void dfs(int u,int s,int w)
{
	if(w>dis[d])
	return ;
	if(u==d)
	{
		if(s>sum&&dis[d]==w)
		{
			res=per;
			sum=s;
		}
		cnt++;
		return ;
	}
	for(int i=0;i<n;i++)
	{
		if(!vis[i])
		{
			vis[i]=1;
			per.push_back(i);
			dfs(i,s+a[i],w+g[u][i]);
			vis[i]=0;
			per.pop_back();
		}
	}
}
int main()
{
	memset(g,0x3f,sizeof g);
	cin>>n>>m>>s>>d;
	for(int i=0;i<n;i++)
	cin>>a[i];
	for(int i=0;i<m;i++)
	{
		int u,v,w;
		cin>>u>>v>>w;
		g[u][v]=g[v][u]=min(g[u][v],w);
	}
	bfs();
	memset(vis,0,sizeof vis);
	per.push_back(s);
	dfs(s,a[s],0);
	cout<<cnt<<' '<<sum<<endl;
	for(int i=0;i<res.size()-1;i++)
		cout<<res[i]<<' ';
	cout<<res[res.size()-1];
	return 0;
} 
